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DESCRIPTION 

APPARATUS AND METHOD FOR NOISE-PREDICTIVE 
MAXIMUM UKELIHOOD DETECTION 

TECHNICAL FIELD 

Present invention relates to data detection apparatus and metliods. It is in particular ap- 
plicable to direct access storage (DASD) devices using partial-response (PR) signaling 
and maxinnum likelihood sequence detection (MLSD). 

BACKGROUND OF THE INVENTION 

In DASD systems as well as in data transmission systems, the signal read from the stor- 
age device or received at the output of a transmission channel has to be converted in the 
receiver into a sequence of symbols which most likely represents the data (symbol) se- 
quence initially stored or transmitted by the sender, despite interference between adja- 
cent symbols, and despite added noise. 

The optimum receiver for detecting an uncoded data sequence in the presence of inter- 
symbol interference (ISI) and noise consists of a whitened-matched filter followed by a 
Viterbi detector which perfonms maximum likelihood sequence detection on the ISI trellis. 
In PRML (Partial-Resj3onse Maximum Likelihood) systems the composite noise at the 
input of the detector is not white, resulting in sub-optimal perfomiance. This sub-op- 
timality is more pronounced at high recording densities since the linear PR4 (PR class 4) 
equalizer tends to enhance the electronic noise component of the composite noise at the 
input of the PRML detector. In addition to colored noise, nonlinear distortion due to non- 
linear bit-shift and head asymmetry could further degrade perfomiance of PRML sys- 
tems. 

Intemational Patent Application PCT/W097/11544 (IBM) discloses a scheme for data 
detection in a direct access storage device, called NPML (Noise-Predictive Maximum U- 
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kelihood detection), that arises by imbedding a noise prediction/whitening filter into the 
branch metric computation of the Viterbi detector. It has been shown via simulations that 
the NPML detectors offer substantial performance gains over the PRML detectors. 

However, Implementation of the Imbedded predictor/whitener as a bank of finite-im- 
pulse-response (FIR) filters or equivalent table look-up operations, e.g. by means of ran- 
dom access memory (RAM), demands a certain level of complexity of the NPML 
hardware realization. A further difficulty with presently available implementation technol- 
ogies is the time constraint imposed by high speed clock requirements. 

Application of the Viterbi algorithm to channels represented by an IIR (Infinite Impulse Re- 
sponse) discrete-time equivalent model has been studied in a publication by A. Duel- 
Hallen et al. "Delayed decision-feedback sequence estimation," IEEE Trans, Commun., 
COM-^7, PP.42&-436, May 1989. The approach in that reference assumes a decision 
feedback equalization prefilter whose feedback filter is an IIR filter combined with Viterbi 
detection. In contrast, the present invention considers a forward partial-response or gen- 
eralized partial-response linear equalizer and a noise whitening IIR predictor imbedded 
into the Viterbi detector. 

OBJECTS OF THE INVENTION 

It is a main object of present invention to provide a more efficient and simpler realization of 
the noise prediction/whitening mechanism for a NPML detector. 

It is a further object of the invention to devise a partial-response maximum likelihood 
detection scheme that utilizes past decisions from the path memory of the Viterbi detector 
for the noise prediction process. 

Another object of the invention is to enable additional mechanisms for a maximum likeli- 
hood detector which allow compensation of reading head asymmetries and of DC offset in 
the receiver. 
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SUMMARY OF THE INVENTION 

The invention provides a partial-response signaling, maximum likelihood detector with 
an infinite impulse response (IIR) noise predictor and whitening filter imbedded in the Vlt- 
erbi detector (INPML). As an additional feature, a DC (zero frequency) notch filter may be 
provided for DC offset compensation. Furthermore, an additional compensation method 
can be provided that is imbedded in the branch metric computation of the INPML detector 
for dynamic signal-asymmetry and DC offset compensation. 

Advantages of the invention are that, besides its simplicity in implementation, it does not 
compromise performance of the detector. It has the further important advantage that it 
can be "piggy-backed" on existing PRML systems so that there is no need for develop- 
ment and implementation of an entirely new channel architecture which would be a com- 
plex and costiy task. The addition of a DC-notch filter or DC compensation method 
renders the INPML detector completely immune to DC offset, tolerant and robust against 
various types of non-linearity and in particular head asymmetry, and tolerant to themnal 
asperities, which is a further impairment in digital magnetic recording systems. 

In the following, embodiments of the invention are described with reference to tiie draw- 
ings. 

LIST OF DRAWINGS 

Fig. 1 is a block diagram of a PRML detector in which present invention can be implement- 
ed. 

Fig. 2A is a block diagram of the maximum likelihood detector according to the invention 
with a single-pole imbedded IIR noise predictor . 

Fig. 2B shows an altemative, equivalent implementation of the IIR noise predictor for the 
inventive maximum likelihood detector of Fig. 2A. 

Fig. 3 shows the 4-^tate trellis for the INPML detector of Fig. 2A. 
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Fig. 4 is a block diagram of an alternative maximum likelihood detector in accordance with 
the invention in which an imbedded IIR noise predictor with 2 poles and 2 zeros is pro- 
vided. 

Fig. 5 shows the bit error probability for various ML detectors including the INPML detec- 
tor of the invention for a channel with the nonnalized linear recording density 
PW5Q/T « 2.86. 

Fig. 6 shows the relative SNR for a given bit error probability, also for various ML detectors 
including INPML (4 states), as a function of PW50/T. 

Fig. 7A is a block diagram of an inventive INPML detector with an additional, imbedded 
DCHfiotch filter. 

Fig. 7B schematically shows a first arrangement of an INPML detector with a cascaded 
DC-notch filter (between equalizer and detector). 

Fig. 7C schematically shows a second arrangement of an INPML detector with a cas- 
caded DC-notch filter (preceding the equalizer). 

Fig. 8 illustrates a RAM based dynamic non-linearity estimation to be used for DC offset 
and head asymmetry compensation. 

Fig. 9 shows the bit en-or probability for the INPML detector with and without asymmetry 
estimation (channel PW50/T - 2.86). 

DETAILED DESCRIPTION 

Fig. 1 is the block diagram of a known PRML channel architecture for a data recording 
system. Such a system was also described in the above mentioned International Patent 
Application PCT/W097/11544. 

Customer data /„ are written in the form of binary digits Un e {-1 , +1 } by write head 15 on 
the disk 11 after being encoded in an encoder 12 by a rate-8/9 RLL (Run Length Limited) 

code and serialized in a serializer 13 and preceded in a precoder 14 by a ^ opera- 
tion where D is the unit delay operator and © means addition modulo 2. When receiving 
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the customer data from said disk 1 1 , an analog signal r(t) is generated by the read head 1 5 
and provided at the read head's output. The signal r(t) is then applied via arm electronic 
circuitry 1 6 to a variable gain amplifier (VGA) circuit 17. The output signal of the the VGA 
circuit 17 Is first low-pass filtered using an analog low--pass filter 18 (LPF) and then con- 
verted to a digital form x„ by an analog-to-digital (A/D) converter 1 9. The functions of the 
A/D converter 19 and VGA amplifier 17 are controlled by the timing recovery and gain 
control loops 20/26 and 21/26, respectively. The analog low-pass filter 1 8 is preferably a 
filter which boosts the higher frequencies to avoid saturation of the A/D converter 19 and 
to support the digital equalizer 22. The digital samples x„ at the output of the fiJD convert- 
er 1 9 (line labeled A in Fig. 1 ) are first shaped to PR4 signal samples by the digital equaliz- 
er 22 (line labeled B in Fig. 1) and are then passed on to a maximum likelihood detector 
circuitry 10 (ML detector) in the form of digital noisy samples yn. The output data of the ML 
detector 1 0 (i.e. the final decisions on the line labeled C in Fig. 1 ), after inverse preceding 
by means of an inverse precoder 23 performing a 1 © operation, are fed via a deser- 
ializer 24 to a rate 9/8 decoder 25 for the rate-8/9 RLL code which delivers the retrieved 

customer data The inverse precoder function 23 following the ML detector in Fig. 1 
can be a separate functional block (as shown) or it can be imbedded In the trellis (survivor 
path memory) of the ML detector. 

In the above mentioned International Patent Application PCT/W097/11544, the maxi- 
mum likelihood detector 10 was a noise-predictive maximum likelihood detector (NPML 
detector) which uses a finite impulse response (FIR) filter/predictor to estimate the noise 
contents of the input signals. 

First Embodiment of invention: 

Present invention provides an improvement for the maximum likelihood detector portion 
of a PRML data recording channel. Its principles are shown in the block diagram of Fig. 2A 
which represents the ML detector portion 10 of Fig. 1. The PR4-equalized samples yn 
appearing on line 30 (marked B) are modified by the noise prediction process according to 
the invention, fed as samples Zn to the Viterbi detector 31 which on its output 32 delivers 
the recovered data symbols a„ which are the delayed final decisions (marked C) of the 
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Viterbi detector. Preliminary symbols (path history decisions) Cn are delivered on an out- 
put 33 of the Viterbi detector. These preliminary symbols are fed to a reconstruction 

means 34 having transfer function 1 - D^, to obtain on line 29 a reconstructed sample 
sequence (a„ - ^n-t)- '^^^^ reconstructed sample sequence is subracted in subtracting 
means 39 from the noisy equalized samples yn appearing on branch line 38, to obtain 
noise samples Wn on line 28. These noise samples Wn are fed to noise predictor circuit 35 
which is an infinite Impulse response (IIR) filter. The predicted noise samples Wn appear- 
ing on output line 36 are combined in circuit 37 with the noisy samples y„ to provide sam- 
ples z„ for input to the Viterbi detector 31. Note that the noise prediction and metric 
update operations are perfonned separately for each state using decisions from the cor- 
responding path memory. These per state operations in the resulting INPML detector are 
symbolized in Fig. 2 by double signal flow lines. 

The IIR predictor 35 of present preferred embodiment is a single-pole device. It receives 
the noise samples vv„ on line 28, which then are delayed in unit delay means 40. Its out- 
put, after combination with a delayed noise signal (output of unit delay means 43) in adder 
41 and multiplication by a coefficient a in multiplication means 42, provides the desired 
noise signal vv^ on line 36. The IIR noise predictor 35 executes the transfer function 
aD/(l - aD). 

Thus, the noise predictor is an infinite impulse response device which is imbedded (for 
each state) into the maximum likelihood detection arrangement, as shown in Fig. 2A. This 
arrangement will be termed INPML detector in the following (for "Infinite impulse response 
Noise Predictive Maximum Likelihood'' detector). Its operation is described in the sequel. 

Let yn be the output of the PR4 digital equalizer, line labeled B in Fig. 1 and Fig. 2A. This 
output then consists of a PR4 data signal, with noiseless nominal values of 
- 2, 0, + 2, and colored noise, i.e., 



where a„ E {-1, + 1} denotes the encoded/precoded data sequence written on the 
magnetic medium with a rate 1/7 and Wn represents the colored noise sequence at the 
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output of the digital partial-response equalizer (T is the channel encoded bit time interval 
and subscript n indicates time instant nT). The power of the colored noise component Wn 
in (1) can be reduced by noise prediction. 

Figure 2A shows the basic concept of a single coefficient MR predictor imbedded Into the 
Viterbi detector. The basic principle can be explained as follows. Let 

^^^^^T^D . «<1 (2) 

denote the D-transform of the single-pole IIR prediction error (whitening) filter. Then, 
assuming ideal path history decisions , i.e., Wn - w„ on line 28, the D-transfonm of 
the whitened noise sequence at the output of £(£>) is given by 

e(D) = 7-z^vv(D) = [1 + Pm w(D) = w(D) + w(D) , 

with 

W = . ,4, 

w(D) in (3) represents the />transform of the colored noise sequence In (1) and w{D) 
denotes its predicted value through the single-pole prediction filter. In the time domain 
w(D) satisfies the following recursive equation 

= a(w^«, + Wn-i) +««-3) = «^n-i • (5) 

The single-pole IIR predictor filter 35 described by (4) and (5) and shown in Fig. 2A can 
also be implemented as shown in Fig. 2B. This IIR predictor 35b also receives noise sam- 
ples Wn on line 28. They are combined in adder 41 with predicted noise samples Wn . then 
delayed by unit delay element (0) 40b and then multiplied in multiplication means 42 by a 
coefficient a. to provide predicted noise samples Wn on output line 36. 

Equation (5) and PR4 signalling define a 4-state IVIL detector which encompasses a 
single-pole IIR predictor. If Sj , y = 0, 1 , 2, 3, denotes the >-th state _ i , _ 2) of the 
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4-state PR4 trellis then, in view of (1) and (3), the branch metric at time nToi the Viterbi 
detector corresponding to the state transition Sj s^^ takes the form 

2 2 

where a„ and a„-2(-s[/) are determined by the hypothesized state transition sj 5^ and 
WniSj) is the predicted noise sample associated with state sj. In particular, the recon- 
structed PR4 samples j - a„_3) in (5) are replaced in practice by (a„_i - 5„,3) 
which are obtained from the path memory of the Viterbi detector associated with state Sj . 

The 4-state INPML detector with a single-pole MR predictor/whitening filter operates on 
the trellis shown in Figure 3. Every state at time instant nT^is associated with : 

-State Metrics: MO), J„(2), MS) , 

- Noise Estimates: iVn(O), w/n(1). w'nCS), Wn{3) . 

The branch metrics shown on the trellis branches of Figure 5 are computed as follows: 
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, where 
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= y„ - 2 + w„(0) 
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, where 
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= yn + H'„(l) 
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, where 


e»(2.3) 


= y« - 2 + w„(2) 
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1 e„(3,3)| 
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, where 


e»(3,3) 


= + w„(3) 



(7) 



The state metrics and noise estimates at the next time instant (n + 1)1 are computed as 
follows: 
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^+i(0) = min { 7„(0) + A„(0.0) . /„(!) + A„(1,0) } 

or 

a e„(l 0) the previous "winning" state 



((\\ - / '"or^ Depending on whether state 0 or state 1 was 



= rnin { J„(2) + A,(2, 1) , y«(3) + A„(3, 1) / 

r ^ ^n{^A) Depending on whether state 2 or state 3 was 

- / %n the previous "winning" state 



y„+l(2) = mm / 7„(0) + A„(0,2) , + A„(l,2) ; 

n\ - f ^ ^^or^^^ Depending on Whether state 0 or State 1 was 

^n^\y^) I ^^^^12) the previous "winning" state 



^n+i(3) = min { U2) + A„(2,3) , 7„(3) + A„(3.3) } 
a e„(2,3) 

n^ = / or Depending on whether state 2 or state 3 was 

W/i+iW / ^ 3j the previous "winning" state 



By setting a = 0 the INPML detector becomes a conventional 4-state PRML detector. 



Second Embodiment of Invention: 



An alternative embodiment of the noise predictor in the inventive INPML detector is now 
described with reference to the block diagram of Fig. 4. In this case the prediction error 
(whitening) filter includes two poles and two zeros. 

The predictor 45 shown in Fig. 4 has the same connections to the other portions of the 
Vlterbi detector as predictor 35 in Fig. 2A: A line 28a (28) transferring noise samples , 
and a line 36a (36) for the predicted noise samples iv„. 

The noise samples w„ appearing on input 28a are twice delayed in unit delay elements 47 
and 48, then multiplied in multiplier 49 by factor [- (Oj + i^j)], and then combined in 
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adder unit 50 with "auxiliaiy signals on lines 51, 52. and 53 to provide predicted noise 
samples )v„ on line 36a. 

The "auxiliary" signals are derived as follows: Signal on line 51 results from multiplying the 
output of delay unit 47 in multiplier 54 by a factor [ - (a , + ^ , )] . Signal on line 52 results 
from multiplying predicted noise samples w,, , delayed in unit delay means 55, by a factor 
[ - a 1 ] in a multiplier 56. Signal on line 53 results from multiplying predicted noise sample 
w„ , twice delayed in unit delay means 55 and 57, by a factor [ - in multiplier 58. MR 
noise predictor 45 executes the transfer function 

Let 



1 + a^D + ajD^ 



EiD) = , (8) 



denote the [>-transform of the two-poles, two-zeros IIR prediction error (whitening) filter. 
Then, the D-^ransfonn of the whitened noise sequence at the output of EiD) is given by 

e{D) = 1 4.^|^^^^2 ^(^) = + ^(^) = ^(^) + ^(^) ' (9) 

with 

^^^^ 1 + a,D + a^D^ ' <^^> 

w(D) in (9) represents the D-transform of the colored noise sequence in (1) and w{D) 
denotes its predicted value through the 2-poles/2-zeros prediction filter. In the time do- 
main vv(D) satisfies the following recursive equation 
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Equation (11) and PR4 signalling can define a 4-state Viterbi detector which encom- 
passes a 2-zeros/2-poles MR predictor. If Sj , / = 0, 1 , 2, 3, denotes the >-th state of the 
4-state PR4 trellis then, in view of (1) and (9), the branch metric at time nToi the Viterbi 
detector corresponding to the state transition Sj f ^ takes the forni 

2 2 

where an and a^^2^sp are determined by the hypothesized state transition sj 5^ and 
Wn{sp Is a noise estimate associated with state sy As in the previous embodiment the 
branch metrics (12) for each state transition can now be computed explicitly from expres- 
sions (7) with WniSj) according to (11). 

Following the same approach, the recursive equation of a 1-zero/2-poles MR predictor 
filter with the corresponding branch metric equation can be readily derived = 0 in 
(8)). Similarly, the recursive equation of a 2-polesllR predictor filter with the correspond- 
ing branch metric equation can also be obtained (^Sj = 0 and = 0 In (8)). 

Performance (Improvement gained by the invention): 

For the magnetic recording channel it has been found that in the entire range of 
0.5 < PW50/T <3.5 where PW50 is the pulse width at the 50% amplitude point of the 
channel's step response and 7 is the duration of the channel encoded bit, II R predictors 
with at most 2 poles and 2 zeros offer the best possible perfomnance irrespective of the 
noise source, i.e., electronic noise or disl< noise, or a combination thereof. In particular, 
for 0.5 ^ PW50/T ^ 2.7 the best possible perfomnance is achieved with an MR predictor 
consisting of 2 poles and 1 zero, whereas for 2.7 < PW50/T < 3.5 best possible perfor- 
mance is achieved with an II R predictor consisting of 2 poles and 2 zeros. If the dominant 
noise source is electronic noise then an INPML detector with a single-pole MR predictor 
offers already substantial perfomnance gains over PRML detectors in a wide range of lin- 
ear recording densities. 
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The error performance of a magnetic recording system employing INPML detection has 
been studied by computer simulation. Figure 5 shows bit error probabilities of simulated 
4-state INPML and NPML detectors for the channel model with PW50/T « 2.86. Curve 
2 shows the perfomnance of INPML4 (4-state INPML) with a single-pole predictor; 
compared to curves 1 (PRML) and 3 (EPRML), a gain in noise margin of 2 dB and 0.5 dB is 
obtained, respectively. Curve 4 corresponds to 4-state NPML with A/ = 4 predictor coeffi- 
cients. Compared to curve 3 an additional gain of 0.4 dB is obtained. Finally, curve 5 
corresponds to a 4-state INPML detector with a 2-zeros and 2-poles IIR predictor. 
Compared to curve 4 an additional 0,2 dB is obtained, and compared to PRML 2.7 dB gain 
is obtained. 

Figure 6 shows the relative SNR (SignaMo-Noise Ratio) required to achieve a bit error 
probability of = 10 " ^ for the 4-state NPML with A/ = 4 predictor coefficients. 4-state 
INPML with a single-pole IIR predictor, and EPRML, all with respect to a PRML system, 
as a function of the nonmalized linear density PWSO/Tfor the Lorentzian channel model. 
Cun/e 1 corresponds to our reference PRML system. Curve 2 con^esponds to the 8-state 
EPRML system and curves 3 and 4 correspond to 4-state INPML with a single-pole ilR 
predictor filter and 4-state NPML/At=4, respectively. We observe that the single-pole 
INPML system is very close in perfomnance to the NPML system and that both provide 
superior perfomnance compared to all the other systems at any linear recording density. 

Further Improvements (Additional Features): 

(a) Including a Notch Filter : 

Further improvement of the invented INPML detector described so far can be achieved by 
incorporating a DC-notch filter for eliminating sensitivity to DC offset. 

The transfer function of a digital notch filter is given by a second order rational function 
with two parameters that characterize the location and width of the notch (see article by P. 
Regalia, S. K. Mitra and P. P. Vaidyanathan The Digital All-Pass RIter: A Versatile Signal 
Processing Building Block," PROCEEDINGS of the IEEE, VOL76. pp. 19-37, January 
1 988). The sensitivity of INPML detectors to DC offset can be eliminated by utilizing such 
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a filter with a notch at DC. It can be seen that the transfer function of a digital DC-notch 
filter is given by 

where 6 is a constant that determines the width of the notch. For a very narrow notch at 
DC, b must be chosen very close to unity. In this case, (1 + b)/2 « 1 and the transfer 
function of the notch filter simplifies to 

^(^) = • (14) 

Digital filters with a notch at DC are very simple to implement and do not require multipli- 
ers. Similarly, the transfer function of a DC-notch analog filter is given by 

ms) = jf-^ . (15) 

where c is a constant that determines the width of the notch. For a very nan^ow notch c 
must be chosen very close to zero. 

The DC-notch filter can be either imbedded into the INPML detector as part of the whiten- 
ing IIR predictor or cascaded with the INPML detector. Figures 7A, 7B and 7C show the 
block diagrams of three possible configurations of an INPML detector with a DC-notch 
filter. In Fig. 7A, the notch filter 44a Is included in the feedback path between the second 
output of the Viterbi detector 31 and its input, the path which also contains the IIR predic- 
tor 35. In the configuration of Fig. 7B, the notch filter 44b is inserted between the digital 
equalizer 22 and the INPML detector 10, and thus receives the samples y„ at its input. In 
the case of Fig. 7C, the notch filter 44c precedes the digital equalizer 22. Clearly, the 
same approach can be applied to other DC sensitive detectors such as NPML, non- 
adaptive DFE (decision feedback equalizer) and FDTS (fixed delay tree search) detec- 
tors. Summarizing, a DC-notch filter renders INPML : 

-completely immune to DC offset, 

-tolerant to MR (magneto-resistive) head asymmetry, 

-tolerant to themnal asperities. 
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(b) Adaptive Asymmetry and DC Estimation Combined with INPML Detection : 

An MR head is a nonlinear transducer that Introduces asymmetry in the readback signal. 
Positive peaks can differ in amplitude from negative peaks by as much as 30 %. DC in- 
sensitive detectors such as PRML and INPML cascaded with a DC-notch filter are quite 
tolerant to small amounts of head asymmetry and suffer only a small perfonmance degra- 
dation. The high losses in perfomiance due to large amounts of head asymmetry can be 
recovered by an asymmetry estimation/compensation algorithm which is imbedded into 
the branch-metric computation of the INPML detector. Consider the following generic 
signal asymmetry model: 



F+2 for = landa„_2 ~ " ^ 
Fq for a„ = a„.2 

^ F^2 for fli, = - landa„^2 = 1 



Then, the branch metric computations of the 4-state INPML detector whose trellis is 
shown in Fig. 3 must be modified as follows: 
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e«(3.3) 


= yn 


- Fo + H>„(3) 



The State metrics and feedback estimates at the next time instant {n + l)rare computed 
as follows: 
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Jn^\(0) = min f MO) + A„(0,0) . y„(l) + A„(1,0) } 

w (0) ^ / ^ ^^or' Depending on whether state 0 or state 1 was 

+ 1 ^ ^ a ( 1 , 0) previous "winning" state 



7„+j(l) = min { U2) + A„(2, 1) , y„(3) + A„(3, 1) } 

a ^„ (2, 1) Depending on whether state 2 or state 3 was 



^ ^ ,a n the previous "winning" State 



y„+,(2) = mm / 7„(0) + A„(0,2) , + A«(l,2) ; 

^ (2) z=z / ^ ^^o?' Depending on whether state 0 or state 1 was 

" + ^ ^ ^ ' a e„ ( 1 , 2) previous "winning" state 



y„+,(3) = min I U2) + A„(2,3) , /„(3) + A„(3,3) ; 

- f 31 = / ^ ^"or^' Depending on whether state 2 or state 3 was 

" + > ^ ^ ' a c„ (3, 3) previous "winning" state 



The values of the nonlinear function /(a„ , a„_2) be estimated for a particular head 
and media combination during manufacturing and then be used as fixed predetermined 
values in the branch metric computation unit of the INPML detector. 

Altemativeiy, an algorithm that estimates these values dynamically can be used. Dynam- 
ic estimation can be achieved by a table look-up stochastic gradient algorithm which up- 
dates one memory location per time interval. Consider delayed decisions from the path 
memory of the INPML detector, say a„.^, fl„.2-rf- These decisions then fomi the ad- 
dress to the memory location of a RAM in which the conresponding values are stored. A 
simple schematic of a RAM with four memory locations containing corresponding values 
at time nTis shown in Figure 8. At the next time instant, (n + 1)7, the content of the 
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memory location addressed by a„_2_^) is updated according to the following 

algorithm: 

^n-d = [ yn-d -//,( K-d^K-l-d) 1 ' 

where y^^^ is the delayed, nonlineariy distorted received signal and is a small 
constant. The above algorithm can also be applied to remedy other forms of nonllnearity. 
In general, the number of symbols which the nonlinear function /(.) depends on deter- 
mines the size of the RAM. 

The following asymmetry model has been used in the simulations: 



2 + 4<5 + DC = 


^+2 


for 


On - landa„_2 = - 1 


0 + DC = 




for 


= «n-2 


2 + 4<J + DC = 


F-2 


for 


fln = ~ landfl„_2 = 1 



This model is basically a second-order nonllnearity coupled with a DC offset. In the simu- 
lations (5 = 0.1, i.e., 33% asymmetry, and DC = 0.05 have been chosen. 

Performance of INPML detector with and without Asymmetry Compensation: 

Figure 9 shows the error perfomnance of a 4-state INPML detector with a single-pole pre- 
dictor in the presence of electronic noise and head asymmetry for the channel model with 
PWSQ/T « 2.86. The head asymmetry is 33%. Curve 1 shows the performance of 
INPML without any asymmetry compensation. Curves 2 and 3 correspond to INPML 
combined with adaptive asymmetry estimation and fi - 0.05 and fi = 0.005, respec- 
tively. Finally, curve 4 shows the perfomnance of INPML with ideal asymmetry estimation/ 
compensation. These performance results show that the dynamic estimation methods 
can achieve the performance of ideal estimates. 
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Digital or Analog implementations: 

Note that, besides a fully digital implementation, the whole or portions of the elements 
shown in the INPML detector 10 of Fig. 2A/B, Fig. 4 and Fig. 7, can be implemented in 
analog circuit technology. ^> 
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CLAIMS 

1 . Detector for encoded and uncoded data sequences having passed a partial-response 
equalizer, including means (10) for maximum likelihood sequence detection with noise 
prediction; 

characterized in that it comprises: 

- an infinite Impulse response (IIR) noise predictor (35; 45) imbedded Into the maximum 
likelihood sequence detection means. 

2. Detector in accordance with claim 1 , in which said infinite impulse response noise pre- 
dictor (35) executes a transfer function aD/{\ - aD), where a is a preselected scaling 
factor and D is the unit delay operator. 

3. Detector In accordance with claim 1 , in which said infinite impulse response noise pre- 
dictor (45) executes a transfer function 

[- (aj + )8,)D - + ^2^D^y^^ + «2^^]' where a,. /Sp and ^2 ^re 

preselected scaling factors and D is the unit delay operator. 

4. Detector in accordance with claim 2, in which 

- said infinite impulse response noise predictor (35) is a single-pole filter comprising: 

- an input (28) for noise samples {w^; an output (36) for predicted noise samples (h>„); 
first unit delay means (40) for delaying the input noise samples (w„); adding means (41) 
for adding the output of the first unit delay means and the output (iv„. {) of a second unit 
delay means (43); and multiplication means (42) for multiplying the output of said adding 
means (41 ) by a constant factor (a), furnishing to the output (36) the predicted noise sam- 
ples {wn)\ the second unit delay means (43) receiving at its input said predicted noise 
samples. 
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5. Detector in accordance with claim 3, in wliich 

- said infinite impulse response noise predictor (45) is a double-pole and double-zero 
filter comprising: 

-an input (28a) for noise samples (vv^ ); an output (36a) for predicted noise samples (vP„); 
first and second unit delay means (47, 48) connected in series, for delaying the input sam- 
ples {^n )* first multiplication means (49) for multiplying the output of said second delay 
means by a first constant factor [ - (02 + fi2^y* adding means (50) for adding the output 
of said first multiplication means and correction terms appearing on further inputs (51 , 52, 
53) of said adding means, furnishing on its output (36a) the predicted noise samples (>v„); 
third and fourth unit delay means (55, 57) connected in series to the output of said adding 
means; and second, third and fourth multiplication means (54, 56. 58) for multiplying the 
outputs of said first (47), third (55) and fourth (57) unit delay means each by a respective 
constant factor ([- (a^ + a,] ;[- a2l)^"^f"r"'s*^'"9^*^®'''0"^P"ts to said fur- 

ther inputs (51, 52, 53) of said adding means (50). 

6. Detector In accordance with claim 4 or claim 5, in which 

- subtracting means (46) are provided for furnishing said noise samples (w„ ) to the input 
(28; 28a) of said infinite impulse response noise predictor (35; 45); 

- said subtracting means (46) having a first input (29) for reconstructed PR4 samples 
(a„ " a^_2) ® second input (38) for equalized, noisy data samples (y„); said sub- 
tracting means (46) forming the difference vv^ of the noisy data samples and the recon- 
structed PR4 samples. 

7. Detector in accordance with claim 6, further including: 

- means (34) for reconstructing PR4 data samples (an - «„-2) "^V ^ transfer function 
( 1 - D^), the input of the reconstructing means being connected to the path history of a 
Viterbi detector (31) providing path history decisions (an )» and its output (29) being con- 
nected to the first input (29) of the subtracting means (46). 
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8. Detector in accordance with claim 1, in which the filter function of said infinite impulse 
response noise predictor (35; 45) has not more than two poles and not more than two 
zeros. 

9. Detector in accordance with claim 1 , further comprising 

- a DC-notch filter (44). either imbedded in, or cascaded with, the maximum likelihood 
detection means (10) which comprises the infinite impulse response noise predictor (35; 
45), for DC offset compensation, 

10. Detector in accordance with one of the preceding claims 1 through 9, in which 

- at least a portion of said maximum likelihood detection means (1 0) which comprises the 
infinite impulse response noise predictor (35; 45), is implemented in analog circuit 
technology. 

11 . Detection method for encoded or uncoded data sequences having passed a partial- 
response equalizer, based on maximum likelihood sequence detection with noise predic- 
tion; 

characterized by: 

- an infinite impulse response (MR) noise prediction procedure integrated with the maxi- 
mum likelihood sequence detection process. 

12. Detection method in accordance with claim 11 , in which said infinite impulse response 
noise prediction filters the difference {w„ ) between partial-response equalized, noisy 
samples {y„) and reconstructed samples (an - -2) derived from the output of a maxi- 
mum likelihood detection process, and provides predicted noise samples (w„). 

13. Detection method in accordance with claim 11 or claim 12, in which 

- the infinite impulse response noise prediction procedure executes a transfer function 
aD/{l - aD), where a is a preselected scaling factor and D is the unit delay operator. 
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14. Detection method in accordance with claim 11 or claim 12, in which 

- the infinite impulse response noise prediction procedure executes a transfer function 
[- (aj + i3,)D - (aj + ^2)^^1/11 + ^jD + ttjl)^], where a^, /Sj, and are 
preselected scaling factors and D is the unit delay operator. 

1 5. Detection method in accordance with claim 1 1 , in which said infinite impulse response 
noise prediction procedure is executed multiple times, once for each state of the trellis of 
the maximum likelihood sequence detection process. 

16. Detection method in accordance with claim 11 , including 

- a procedure for DC offset compensation and reading head asymmetry compensation, 
comprising following steps: 

A 

- addressing a memory containing estimates of a nonlinearity, using delayed deci- 
sions a^_^ as addresses, and using the nonlinearity signal estimate values thus ob- 
tained from the memory, in the branch metric computation for each state of the maximum 
likelihood sequence detection process which includes the infinite impulse response noise 
prediction procedure; and 

- concurrently, updating the estimates of the nonlinearity j^.) contained in the memory, 
by a stochastic gradient procedure according to 
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